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ABSTRACT 

We present Chandra snapshot observations for a sample of 7 sources selected from the 
Arecibo OH megamaser (OHM) survey at z ~ 0.13 — 0.22 and with far-infrared luminosities 
in excess of 10 11 L©. In contrast with the known H2O megamasers, which are mostly asso- 
ciated with powerful Active Galactic Nuclei (AGN), the situation is far less clear for OHMs, 
which have been poorly studied in the X-ray band thus far. All of the observed sources are 
X-ray weak, with only one OHM, IRAS FSC 03521+0028 (z = 0.15), being detected by 
Chandra (with 5 counts). The results from this pilot program indicate that the X-ray emis- 
sion, with luminosities of less than w 10 42 erg s , is consistent with that from star formation 
(as also suggested in some cases by the optical spectra) and low-luminosity AGN emission. 
If an AGN is present, its contribution to the broad-band emission of OHM galaxies is likely 
modest. Under reasonable assumptions about the intrinsic X-ray spectral shape, the observed 
count distribution from stacking analysis suggests absorption of f 
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Extragalactic hydroxyl (OH) megamaser (OHM) emission has been 
studied since the early eighties when it was discovered in Arp 220 
(Baan, Wood & Haschick 1982). To be produced, OHM activ- 
ity requires (1) high molecular density (n#, = 10 4 ~ 7 cm -3 ; Baan 
1991), (2) a "pump" to invert the hyperfine population of the OH 
ground state, and (3) a source of 18 cm continuum emission to 
stimulate maser emission (Burdyuzha & Romberg 1990), whose 
main lines are at 1667 MHz and 1665 MHz and have luminosi- 
ties Loh=10 1-4 Lq. The galaxy merger environment is able to 
supply all of these requirements: the merger/interaction concen- 
trates molecular gas in the merger nuclei, creates strong far-infrared 
(FIR; 8-1000 /jm) dust emission from reprocessed starburst and ac- 
tive galactic nucleus (AGN) activity, and produces radio continuum 
emission. The FIR radiation field and/or collisional shocks in the 
molecular gas can invert the OH population and thus allow maser 
emission. 

Not surprisingly, all known OHMs have been observed in lu- 
minous infrared galaxies (LIRGs; Lfir> 10 L@), often favouring 
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the most FIR luminous (e.g., Baan 1991; Baan, Salzer & LeWinter 
1998; Darling & Giovanelli 2002a, hereafter DG02a; Lo 2005), the 
ultraluminous infrared galaxies (ULIRGs; Lfir> 10 12 Lq). The 
high FIR luminosities of ULIRGs are commonly thought to arise 
from dust absorption and FIR re-emission of an intense but ob- 
scured starburst and/or AGN radiation field. While it was known 
early on that ULIRGs nearly always show evidence for interactions 
(i.e., collisions/mergers; e.g., Sanders & Mirabel 1996; Clements 
et al. 1996; Borne et al. 2000; Farrah et al. 2001, 2003), many in- 
vestigators have debated the luminosity dependence of the fraction 
of interacting systems among the ULIRG population and the possi- 
bility that these sources may represent the transition stage between 
galaxy mergers and quasars (e.g., Farrah et al. 2001; Tacconi et al. 
2002; see also Yun et al. 2004). 

OHM studies provide powerful diagnostics of the physical 
conditions in the innermost regions of luminous FIR galaxies. 
There is evidence that the OH emission is often produced by an en- 
semble of many masing regions in the nuclear regions of (U)LIRGs 
on scales of a few hundred parsecs or less (Diamond et al. 1999). 
VLBI observations of a few nearby OHMs have demonstrated that 
OH maser emission can arise in circumnuclear discs or tori such as 
for III Zw 35 (Pihlstrom et al. 2001) or Mrk 231 (Klockner, Baan 
& Garrett 2003), or, as for Arp 220 (Rovilos et al. 2003), the emis- 
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Table 1. Chandra Cycles 4 and 5 observations of OH megamaser galaxies: observation log. 



Src. Name 






Optical 


X-ray Obs. 


Exp. Time 


Reference 


IRAS FSC 


z 




«2000 


§2000 


Date 


(ks) 


(Arecibo Survey) 


01562+2528 


0.1658 


01 


59 02.61 


+25 42 35.4 


2003 Nov. 08-09 


4.00 


(1) 


02524+2046 


0.1815 


02 


55 17.09 


+20 58 56.5 


2003 Nov. 08 


4.31 


(1) 


03521+0028 


0.1522 


03 


54 42.21 


+00 37 03.8 


2002 Dec. 25 


9.78 


(1) 


08201+2801 


0.1680 


08 


23 12.62 


+27 51 39.6 


2003 Nov. 20 


4.41 


(2) 


08279+0956 


0.2085 


08 


30 39.38 


+09 46 36.1 


2004 Jan. 13 


4.38 


(2) 


09531 + 1430 


0.2151 


09 


55 50.11 


+ 14 16 07.9 


2003 Dec. 09 


4.03 


(2) 


09539+0857 


0.1290 


09 


56 34.30 


+08 43 06.1 


2004 Jan. 06 


4.81 


(2) 



NOTES — The optical positions have been derived from the Digital Palomar Sky Survey (DPOSS2) fi-band images 
using SEXTRACTOR (Bertin & Amouts 1996). REFERENCES — (1) Darling & Giovanelli 2002a (DG02a); (2) 
Darling & Giovanelli 2001 (DG01). 



sion can show a complicated, irregular morphology and velocity 
structure. Blueshifted/redshifted line components in several OHMs 
(Baan, Haschick & Henkel 1989; Pihlstrom et al. 2005) have been 
interpreted as due to bulk outflows/inflows of molecular gas; this is 
expected given the tidal streaming and starburst-driven winds asso- 
ciated with major mergers. It is noteworthy that most VLBI studies 
of OHMs show unresolved maser spots of surprising velocity width 
(tens of km s~'), often attributed to turbulence and suggestive that 
compact OHM emission may not be pumped by the FIR radiation 
field alone (e.g. Lonsdale et al. 1998). A notable exception is the 
infrared quasar Mrk 231, which lacks a compact maser component 
contrary to expectation (Lonsdale et al. 2003; Klockner, Baan & 
Garrett 2003). The relationship of OHMs to starbursts versus AGN 
thus remains murky: all manner of masing is seen coupled with 
all manner of optical classifications. To obtain deeper insight into 
the OHM phenomenon, we require a minimal-obscuration means 
to quantify the role of AGN versus starburst in masing regions. The 
radio-FIR relation hints at some difference between masing and 
non-masing (U)LIRGs (DG02a), but there is a degeneracy between 
starbursts and AGN in this relation. 

In contrast with H2O megamasers, which are likely pumped 
by an AGN X-ray radiation field (e.g., Baan 1997; Braatz, Wil- 
son & Henkel 1997; Henkel et al. 1998, 2005; Townsend et al. 
2001; Maloney 2002; Braatz et al. 2003, 2004), the AGN content of 
OHMs has not been investigated with a minimal-obscuration probe 
of AGN activity, aside from radio continuum observations of radio- 
loud AGN that lie well off of the radio-FIR relation. A previous op- 
tical study of OH maser sources (most of which were OHMs) found 
comparable fractions of AGN- and starburst-dominated galaxies 
(Baan et al. 1998) and a non-negligible fraction of composite spec- 
tra (i.e., showing evidence of both AGN and starburst activity) 
in the OH maser population. Because of their penetrating nature, 
X-rays represent an efficient tool to provide constraints on the en- 
gine powering OHMs. Furthermore, X-rays typically provide maxi- 
mal contrast between the AGN emission and that of the host galaxy 
and/or starburst component (e.g., Vignati et al. 1999). 
Hereafter we adopt Hq=15 km s _1 Mpc~' in a A-cosmology 
(n M =0.3 and Q. A =0J, Spergel et al. 2003). 



2 SAMPLE SELECTION AND SOURCE PROPERTIES 

The goal of the present paper is to provide basic constraints on 
the X-ray emission of OHM galaxies using Chandra snapshot ob- 
servations. Past studies have shown evidence for a correlation be- 



tween the OH line width and the X-ray luminosity in the soft band 
(Kandalyan 2003) using a compilation of OH masers from the lit- 
erature (mostly from Baan et al. 1998) detected by ROSAT. This 
was interpreted as suggesting that X-ray heating of molecular gas 
provides additional collisional excitation to the maser emission, 
although further studies are required to investigate this issue. We 
note, however, that the compilation of OHMs chosen by Kandalyan 
(2003) is heterogeneous, since it was not drawn from a well defined 
sample. For these reasons, we prefer to focus, in this pilot Chan- 
dra program, on the sample of OHMs selected from the Arecibo 
Observatory OH megamaser survey which is a flux-limited survey 
designed to quantify the relationships between merging galaxies 
and the OHMs that they produce, with the goal of using OHMs 
as luminous tracers of mergers (and, possibly, of star formation; 
Townsend et al. 2001) at high redshifts (Darling & Giovanelli 2000, 
2001, 2002b, hereafter DG00, DG01, and DG02b, respectively). 
This survey has been conducted over about one-quarter of the sky 
to a distance limit of roughly 1 Gpc. Candidates for observations 
with the Arecibo Observatory were selected from the PSCz (Saun- 
ders et al. 2000), which is a flux-limited (IRAS f 60 ^ m > 0.6 Jy) 
redshift survey of ~ 15,000 IRAS galaxies over 84 per cent of the 
sky. The Arecibo survey detected 50 new OHMs in the redshift 
range 0.11-0.27 (DG02a,b), thus doubling the sample of known 
OHMs and increasing the z > 0. 1 sample sevenfold. Thanks to this 
survey, it is now evident that the fraction of OHMs in (U)LIRGs 
is an increasing function of Lfir : about one-third of the "warmest" 
(i.e., with higher dust temperatures) ULIRGs are characterized by 
OH megamasers (DG02a). 

The OHM galaxies proposed for observation with Chandra 
have been chosen using the catalog published by DG01 available at 
the time of Chandra Cycle 5 proposals (6 of the 7 OHMs described 
in the following); our targets are among the highest luminosity 
members (Loh > 10 3 L©; ~ 30-40 per cent of the OHM popu- 
lation in the Arecibo survey) of the maser population (filled circles 
in Fig. 1). They also constitute a complete sample, since they were 
targeted in order of right ascension. One additional OHM, charac- 
terized by lower OH luminosity (filled triangle in Fig. 1), has been 
included using an archival Chandra observation. It is worth not- 
ing that our OHMs lie at higher redshifts than most of the OHMs 
with X-ray detections in the literature (e.g., Arp 220 at z = 0.018, 
Mrk 231 at z = 0.042, etc.). On the basis of the merger/interaction 
scenario for OHMs discussed in §1, the present sample, with all 
of the 7 OHMs being characterized by double nuclei 1 and/or tidal 

' We note that, given the presence of double nuclei, the optical positions 
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Figure 1. Logarithm of the OH luminosity vs. FIR luminosity (both in units 
of Lq) for the Arecibo OHM final sample (DG02b). The filled symbols indi- 
cate the OHMs with X-ray observations, while the filled triangle shows the 
X-ray detected OHM. The short-dashed lines indicate the permitted range 
of Lpir when the flux density at 100 /.im is only available as an upper limit 
(see DG02b); in the case of OHMs with X-ray observations, both circles 
are filled. Note that all of the OHM galaxies have Lfir in the LIRG regime, 
with a significant fraction (rs 35 per cent) being ULIRGs. 

tails (Darling et al., in preparation), can be considered representa- 
tive of the OH megamaser population overall. Furthermore, their IR 
colour (log /100//60) distribution shows no significant evidence for 
any observational bias toward "warm" (i.e., lower /100//6O values, 
possibly related to a contribution from an active nucleus) or "cold" 
galaxies (see Fig. 2). 

2.1 Notes on individual sources and optical spectra 

In the following, we report on the OH line and optical properties 
of the 7 OHM galaxies with Chandra observations. A detailed de- 
scription of their optical spectra is presented in Darling et al. (in 
preparation). 

01 562+2528: this source shows the presence of two nuclei (one 
much brighter than the other) in the optical (with a separation of 
6.3", i.e., 17.9 kpc), suggestive of a multiple merger. The East- 
ern component has a starburst classification, while for the Western 
component there is not enough signal for a spectral classification 
(see Table 2 for references). 

02524+2046: this object has the most unusual spectrum of the 
Arecibo OHM survey, showing multiple strong narrow components 
in both OH lines. The optical morphology is elliptical-like with a 

reported in Table 1 are intermediate between the two nuclei; in these cases, 
both nuclei fall in the same 2"-radius aperture used for X-ray manual pho- 
tometry (see §3). For IRAS FSC 01562+2528, we report the position of the 
optically most luminous nucleus, given the association of the OH and radio 
emissions with this nucleus. 



Figure 2. Logarithm of the OH luminosity (in units of Lq) vs. IR colour 
(log /100//6O) for the OHM galaxies from DG02b. Symbols are the same 
as in Fig. 1 . 

single tidal tail (Darling et al., in preparation), while the optical 
spectrum is typical of a starburst galaxy. 

03521+0028: this source is comprised of a close pair; both 
components, separated by 1.3" (corresponding to 3.4 kpc), are visi- 
ble in the Sloan Digital Sky Survey (SDSS) but not in the DPOSS2. 
This object has been classified as a starburst according to the ISO 
spectrum (Lutz, Veilleux & Genzel 1999) and a low-ionisation nu- 
clear emission-line region galaxy (LINER) according to the optical 
spectrum (Nagar et al. 2003). 

08201 +2801 : the HST image shows two nuclei (with a separa- 
tion of 1.1", corresponding to 3.2 kpc), having a starburst classifi- 
cation, and a tidal tail, suggesting an advanced merger (DG01). The 
NVSS flux density reported in DG01 is associated with a source at a 
distance of ~ 23". Nagar et al. (2003) reported a flux-density upper 
limit of 0.8 mJy at 15 GHz, roughly corresponding to a 1.4—15 GHz 
spectral slope of ^ 0.75. 

08279+0956: the optical image of this source (classified as a 
LINER) shows the presence of 3 tidal tails (Darling et al., in prepa- 
ration). 

09531 +1 430: this galaxy is an interacting system whose com- 
ponents are blended in the images from the Automatic Plate Mea- 
suring (APM) facility; the 2MASS position is on the Western com- 
ponent (having a Seyfert 2 spectrum), while the FIRST position is 
on the Eastern component (characterized by a starburst spectrum; 
Darling et al., in preparation). 

09539+0857: the optical image shows the presence of tidal 
tails and a single nucleus, classified as a LINER (Veilleux, Kim 
& Sanders 1999) or Seyfert 2 (Darling et al., in preparation). The 
two classifications were performed in an identical manner. In both 
cases, the [O III]/HP ratio lies very close to the (somewhat arbi- 
trary) dividing line between LINERs and Seyfert 2s, and the mea- 
sured line ratios obtained by both groups agree to within the mea- 
surement uncertainties. Nagar et al. (2003) reported a flux density 
of 1.2 mJy at 15 GHz, corresponding to a 1.4-15 GHz spectral 
slope of ks 0.6-0.8. 



Table 2. Multi-wavelength properties of the Arecibo sample of OHMs observed by Chandra. 



Source 


z 


«H 


R 


fi 


/12 


/25 


/60 


A oo 


logi-FIR 


logiOH 


fl .4 GHz 


/l.4GHz 


F FB 


logZ.2-10 


Class. 


Ref. 










APM 




IRAS 










FIRST 


NVSS 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


01562+2528 


0.1658 


7.46 


15.58 


17.23 


< 0.126 


< 0.133 


0.809 


1.623 


11.90 


3.25 




5.9 


< 11.3 


<4.4x 10 41 


STB 


(1) 


02524+2046 


0.1815 


11.3 


18.08 


20.02 


< 0.055 


< 0.079 


0.957 


< 4.789 


11.81-12.28 


3.74 




2.6 


< 11.4 


< 5.4 x 10 41 


STB 


(1) 


03521+0028t 


0.1522 


12.6 


18.44 


> 21.13 


< 0.250 


0.233 


2.638 


3.833 


12.28 


2.44 




6.1 


3 9 +2.7 


1.3 xlO 41 


STB/L 


(1,2,3,4) 


08201+2801t 


0.1680 


3.59 


17.12 


18.39 


< 0.085 


< 0.162 


1.171 


1.429 


12.00 


3.45 


4.74 




<4.4 


< 1.8 x 10 41 


STB 


(4,5) 


08279+0956 


0.2085 


3.80 


18.08 


20.04 


< 0.079 


< 0.192 


0.586 


< 1.263 


11.75-12.01 


3.23 


2.84 


3.7 


<9.5 


<6.1xl0 41 


L 


(II 


09531+1430 


0.2151 


3.07 


17.25 


> 19.27 


< 0.140 


< 0.169 


< 0.776 


1.040 


12.08 


3.42 


3.23 


2.5 


< 10.2 


< 1.7 x 10 42 


S2/Ltt 


(II 


09539+0857t 


0.1290 


3.08 


17.74 


19.62 


<0.154 


< 0.147 


1.438 


1.044 


11.79 


3.48 


5.32 


8.50 


<8.5 


< 1.9 x 10 41 


S2/L 


(1.3) 



Luminosities are computed using Hq=75 km s Mpc , O. M =0.3 and O. A =0.7. 
t Also in the SDSS Data Release 3 catalog (DR3; Abazajian et al. 2005). 

ft The Eastern nucleus is classified as a starburst, while the Western nucleus is classified as a Seyfert 2 (see text for details). 

(1 ) IRAS FSC Name; (2) redshift; (3) Galactic column density from Dickey & Lockman (1990), in units of 10 20 cm" 2 ; (4-5) R- and fi-band magnitudes from APM; (6-9) IRAS flux densities, in units of Jy; 
(1 0) log of the FIR luminosity (in units of L©) computed according to LpiR=3.96 X 10 5 X D£ (2.58 X / 60 + /loo), where £>l is the luminosity distance in Mpc (see Fullmer & Lonsdale 1989); when /ioo is an 
upper limit, the permitted range of Lfir is reported (see §2.2 of DG02b); (11) log of the OH luminosity, in units of L©; (12) integrated 1.4 GHz flux density from FIRST (Becker, White & Helfand 1995); (13) 
1.4 GHz flux density from NVSS (Condon et al. 1998); (14) Galactic absorption-corrected flux in the observed 0.5-8 keV band, in units of 10~ 15 erg cm -2 s -1 ; (15) 2-10 keV rest-frame luminosity corrected 
for the effects of Galactic absorption and computed from the full-band flux (or upper limit) assuming F = 2.0, in units of erg s~' ; (16) spectroscopic classification: STB=starburst; L=LINER; S2=Seyfert 2; (17) 
reference for the spectroscopic classification. 

List of references — (1) Darling et al., in preparation; (2) Lutz, Veilleux & Genzel 1999; (3) Veilleux, Kim & Sanders 1999; (4) Nagar et al. 2003; (5) Kim, Veilleux & Sanders 1998. 



3 X-RAY OBSERVATIONS AND DATA REDUCTION 

Six out of the 7 OHMs presented here were targeted by Chandra 
during Cycle 5. IRAS FSC 03521+0028, found at the aim point 
of a Chandra archival observation, was observed in Cycle 4. This 
source was targeted in a Chandra program aimed at observing a 
sample of ULIRGs (see Teng et al. 2005); we are confident that the 
inclusion of this object does not bias our "original" sample in terms 
of AGN content (see §4). The observation log is shown in Table 1. 

All of the sources were observed with the Advanced CCD 
Imaging Spectrometer (ACIS; Garmire et al. 2003) with the S3 
CCD at the aimpoint. Standard data reduction was adopted (see §2 
of Vignali et al. 2005 for a detailed description) using the Chandra 
Interactive Analysis of Observations (CIAO) Version 3.2 software. 
Source detection was carried out with WAVDETECT (Freeman et 
al. 2002) using a false-positive probability threshold of 10 -6 . We 
searched for Chandra sources within fa 0.8" of the optical position. 
IRAS FSC 03521+0028 is the only detected source (filled triangle 
in Fig. 1), with 5 counts in the observed 0.5-8 keV band; we note 
that this is also the source with the highest exposure by Chandra 
(see Table 1). Given the adopted threshold and the small number 
of pixels being searched due to the known source positions and the 
sub-arcsec on-axis angular resolution of Chandra, the probability 
that this detection is spurious is extremely low (~ 8 x 10~ 6 ). All of 
the remaining sources (shown as filled circles in Fig. 1) were not 
detected also using higher (i.e., less conservative) threshold values 
(up to 10- 4 ). 

Table 3 summarizes the X-ray photometric results in the soft 
band (0.5-2 keV), the hard band (2-8 keV), and the full band 
(0.5-8 keV) using circular apertures centred on the optical source 
position and an extraction radius of 2". Although two X-ray pho- 
tons falling in adjacent pixels can be considered in some circum- 
stances a detection (given the low background level of Chandra 
observations; see, e.g., Vignali et al. 2005), we prefer to follow 
a more conservative approach and report these sources (2 in the 
present sample; see Table 3) as X-ray undetected. 

The final catalog of OHMs from the Arecibo survey (DG02b) 
was also cross correlated with the ROSAT and XMM-Newton 
archives; we obtained no further X-ray detections. X-ray con- 
straints from ROSAT observations (available for 9 OHMs) are gen- 
erally loose because of the large off-axis angles of the sources 



Table 3. X-ray counts and full-band count rates. 







X-ray Counts 




CR 


Source 


(0.5-2 keV) 


(2-8 keV) 


(0.5-8 keV) 


(0.5-8 keV) 


01562+2528 


<3.0 


<6.4 


< 6.4* 


< 1.60 


02524+2046 


<4.8 


<4.8 


<6.4 


< 1.49 


03521+0028 


3 0+ 2 - 9 

^•"-1.6 


<6.4 


4 0+34 


50+ 35 


08201+2801 


<3.() 


<3.0 


<3.0 


<0.68 


08279+0956 


<6.4 


<3.0 


<6.4* 


< 1.46 


09531 + 1430 


<3.0 


<6.4 


<6.4 


< 1.59 


09539+0857 


<6.4 


<3.0 


<6.4 


< 1.33 



* The two X-ray counts are contiguous. 

Errors on the X-ray counts were computed according to Gehrels (1986). The 
upper limits are at the 95% confidence level and were computed according 
to Kraft, Burrows & Nousek (1991). For the sake of clarity, upper limits 
of 3.0, 4.8, and 6.4 indicate that 0, 1, and 2 X-ray counts, respectively, 
have been found within an extraction region of radius 2" centred on the 
position of the OHM source (considering the background within this source 
extraction region to be negligible). Count rates (CR) are in units of 10~ 3 
counts s _1 . 

and the shallow exposures of the pointed PSPC and HRI obser- 
vations; hence these constraints will not be used in this paper. One 
OHM was targeted by XMM-Newton (IRAS FSC 13218+0552 at 
z = 0.205, with net exposures of fa 8.6-9.7 ks in the EPIC MOS/pn 
cameras) but only an upper limit, consistent with the values we ob- 
tained for our sources (see §4 and Table 2), was derived (Bianchi et 
al. 2005). Recently, it has been suggested that this source, charac- 
terized by an extremely red infrared continuum (Low et al. 1989) 
and a blueshifted component in the [O III] line (e.g., Zheng et 
al. 2002), might harbour a Compton-thick AGN (i.e., with Nu > 
10 24 cm -2 ), surrounded by a strong starburst (see §4.4 and Fig. 6 
of Bianchi et al. 2005). We also note that IRAS FSC 13218+0552 
has the warmest /100//6O colour in Fig. 2 (i.e., it is the left-most 
object in the plot), suggestive of an AGN contribution at IR wave- 
lengths. 




JQ41 IQK 1Q 43 1Q 44 J Q45 J Q46 1 47 

0.5 — 8.0 keV luminosity (erg s _1 ) 



Figure 3. Rest-frame FIR luminosity vs. X-ray luminosity for the OHM 
galaxies presented in this paper (filled circles: X-ray upper limits; filled tri- 
angle: X-ray detected OHM; filled square: stacked X-ray counts from the 
whole sample of 7 OHMs observed by Chandra, see §4.1). The open tri- 
angle indicates the OHM from the Arecibo sample which was observed 
but not detected by XMM-Newton (see Bianchi et al. 2005 and §3). For 
comparison, some other sources are reported (after correcting, when pos- 
sible, for intrinsic absorption in the X-ray band): the sub-mm galaxies in 
the Chandra Deep Field-North (Alexander et al. 2003, 2005a,b) likely con- 
taining an AGN (open circles) and a "pure" starburst (small open squares); 
well known AGN-dominated (mostly in the right-hand side of the plot) 
and star-formation dominated sources (in the left-most region of the fig- 
ure) from the literature plus two distant well-studied sub-mm galaxies 
(SMM 02399-0134 and SMM 02399-0136), all shown as large open 
squares with labels. The slanted lines show ratios of constant X-ray to 
FIR luminosity; this ratio is typically w 3-30 per cent for AGN-dominated 
sources. Adapted from Alexander et al. (2003); see Alexander et al. (2005a) 
for details. Note that some of the sources from the literature sample are low- 
redshift OHM galaxies (marked by a "*" after the source name). 

4 WHAT ARE X-RAYS TELLING US ABOUT THE 
NATURE OF THE OH MEGAMASERS? 

The broad-band properties of the OHM sample targeted by Chan- 
dra are reported in Table 2. Both the photon index and X-ray lu- 
minosity can provide constraints on the nature of the high-energy 
emission from OHM galaxies. However, such constraints, with the 
current data, are basic due to the low counting statistics. All of the 
X-ray sources have rest-frame 2-10 keV luminosities at most of 
~ 10 42 erg s _1 (see Table 2), which is consistent with both star for- 
mation (e.g., Zezas, Alonso-Herrero & Ward 2001; Alexander et al. 
2002) and low-luminosity AGN (LLAGN) activity (e.g., Terashima 
et al. 2002; Terashima & Wilson 2003). The possibility that the ob- 
served X-ray emission is largely due to star formation finds support 
from Fig. 3, where the FIR luminosity is plotted vs. the 0.5-8 keV 
luminosity. Our sources are located in the left-most part of the dia- 
gram, where either starbursts (e.g., Arp 220) or starburst-dominated 
(i.e., with little apparent contribution from the AGN to the bolo- 
metric emission; see Alexander et al. 2005a) sub-mm sources typ- 
ically lie. However, some X-ray obscured AGN whose broad-band 
emission is dominated by star formation are also present in this re- 



gion; a fraction of these sources are lower redshift OHM sources 
(e.g., Mrk 273 and UGC 5101); the possibility that our sources 
could contain absorbed AGN will be discussed in detail using X-ray 
stacking analysis (see §4.1). 

Optical/near-infrared (NIR) spectroscopy (see §2.1) indicates 
that 50 per cent of our sources show the typical signatures of star- 
burst galaxies, 25 per cent have a LINER classification, and 25 per 
cent a Seyfert 2 classification (according to the classification re- 
ported in Darling et al., in preparation; see also column (16) in Ta- 
ble 2); IRAS FSC 09531 + 1430 is counted twice because both nu- 
clei were spectroscopically identified (see Darling et al., in prepa- 
ration). We note, however, that the optical classification of these 
sources is difficult: when more than one nucleus is present, un- 
less the projected physical separation of the nuclei is larger than 
the resolution of the spectrograph, optical spectra will be a blend 
of nuclei, dominated by the optically dominant (i.e., the least ob- 
scured and/or most optically luminous in emission lines) nucleus. 
However, OHMs seem to favour the dustiest environments and may 
select the optically "subordinate" nucleus. Therefore, there may be 
little correlation between OHM properties and optical spectral type. 
Furthermore, we note that optical and NIR studies are sometimes 
not adequate to unveil the engine of the X-ray emission. An em- 
blematic, though perhaps extreme, case of different classification 
at different wavelengths is provided by the ULIRG NGC 6240, 
which is classified as a starburst in the NIR (Genzel et al. 1998), 
a LINER in the optical (Veilleux et al. 1995) and a Compton-thick 
AGN in the hard X-ray band (Vignati et al. 1999; also see Komossa 
et al. 2003). Some further cases similar to that of NGC 6240 are 
presented by Maiolino et al. (2003). 

Finally, it is interesting to note that the X-ray luminosity and 
the X-ray-to-FIR luminosity ratio of the only X-ray detected OHM 
in the present sample do not appear to be different from those of 
the other sources; given the number of observed counts (see Ta- 
ble 3), the detection of this source has been possible because of the 
exposure time which is double that of the other OHMs. 

4.1 X-ray stacking analysis 

We can place tighter constraints on the average X-ray properties 
of the OHMs under investigation by "stacking" their counts. Our 
sources have St^'g counts in the 0.5-2 keV band and counts 
in the 2-8 keV band. Monte Carlo simulations in regions close to 
our targets (see Vignali et al. 2005 for details about the adopted 
procedure) indicate that the background in all of the 7 fields does 
not contribute more than 0.14 and 0.31 counts in the soft and 
hard bands, respectively; therefore we are confident we obtain se- 
cure detections in both bands by stacking the counts from the 7 
OHMs. 2 The ratio of the observed hard vs. soft-band source counts 
(0.85+q 43) suggests that our sources have, on average, a photon in- 
dex flatter (~ 0.7+q g) than 1=1. 8-2.1, which is typically used for 
AGN and starburst galaxies at hard X-ray energies (e.g., Piconcelli 
et al. 2005; Ptak et al. 1999). 3 If the hardness of our sources is due 
to the presence of intrinsic absorption, the implied average column 
density is ~ 1.5 x 10 22 cm~ 2 (ranging between ~ 5.8 x 10 21 cm~ 2 

2 The Poisson probability of obtaining 8 (7) counts or more when 0.14 
(0.31) counts are expected is w 3.2 x 10~ 12 4.2 x 10~ 8 ). 

3 We note, however, that a large population of neutron star high-mass X-ray 
binaries, whose presence is expected in cases of galaxies undergoing intense 
star formation, can produce a flat (r Ri 0.5-1) photon index; see Colbert et 
al. (2004). 



and fa 2.6xl0 22 cm -2 ) at the average redshift of z = 0.17. The 
corresponding absorption-corrected rest-frame 2-10 keV luminos- 
ity would be in the range fa 1-3 x 10 41 erg s~ 1 . Intrinsically flatter 
X-ray slopes would produce lower column densities (by ~ 20 per 
cent, if r = 1.6 instead of 2.0 is adopted), while the presence of an 
additional soft X-ray component (e.g., thermal emission from the 
starburst) would produce the opposite effect. 

Figure 3 shows that, if an active nucleus is present in the 
OHMs under investigation, its contribution to their broad-band 
emission is probably modest, unless the AGN is Compton-thick. In 
the Compton-thick case, the X-ray emission in the Chandra band 
would be the small (presumably a few per cent; see, e.g., Comastri 
2004; Guainazzi, Matt & Perola 2005) scattered/reflected fraction 
of the nuclear X-ray emission, which would be ^ a few x 10 43 erg 
s" 1 (i.e., at most in the Seyfert regime). The reprocessing of the 
X-ray and ultraviolet emissions at longer wavelengths would still 
be insufficient to produce the observed FIR luminosity without 
requiring a dominant contribution from star formation. We note 
that some sources from the literature sample are characterized by 
Compton-thick absorption (e.g., Mrk 231, Gallagher et al. 2002, 
2005, Braito et al. 2004; UGC 5101, Maiolino et al. 2003). 

Although the presence of obscuration is uncertain and needs to 
be confirmed by deeper observations with Chandra for an extended 
sample of OHMs and, possibly, with XMM-Newton for the most 
"promising" sources, the column density allowed by the current 
data is consistent with the X-ray absorption observed in many H2O 
megamasers. These sources, however, are mostly associated with 
X-ray bright (fa 10 42 " 43 erg s" 1 ) and obscured (/% > 10 22 cm" 2 ) 
Seyfert 2 galaxies (e.g., Maloney 2002; Braatz et al. 2003), where 
the amplifying water vapor molecules are most likely pumped by 
collisional processes resulting from X-ray heating in parsec-scale 
molecular tori with densities of fa 10 7 " 8 cm" 3 (e.g., Neufeld, 
Maloney & Conger 1994). Column densities of fa 10 22 cm" 2 are 
not uncommon in local (U)LIRGs (e.g., Risaliti et al. 2000) and 
LLAGN (e.g., Terashima et al. 2002; Terashima & Wilson 2003). 
Given the presence of double nuclei and/or indications of past 
mergers in the OHMs of our sample, it is possible that the starburst 
itself can produce the obscuration, as suggested by some authors 
(e.g., Fabian et al. 1998; Ohsuga & Unemura 2001; Wada & Nor- 
man 2002). 

4.2 X-ray observations of "bona-fide" OHMs from the 
literature 

We obtained further constraints on the X-ray emission from OHM 
galaxies using the sample of 54 "bona-fide" OH masers [both kilo- 
(Loh^ 10 1 Lq) and mega-] from the literature (mostly from Baan 
et al. 1998), listed in Tables 7 and 8 of DG02a. We caution the 
reader against over-interpretation of the results we present in this 
section: this sample is not well defined, lacks either published OH 
spectra or measurements of the 1667 MHz OH line in about 50 
per cent of the cases, and includes some suspect OH detections, as 
pointed out by DG02a. 

Twenty-four OH masers have accessible pointed X-ray obser- 
vations in the ROSAT, Chandra, and XMM-Newton archives (for 
one additional source, the XMM-Newton data are not public yet); 
all of these galaxies are OHMs. Seventeen sources are X-ray de- 
tected; the high detection rate (fa 70 per cent) is due to the fact that 
most of these sources were targeted by the sensitive X-ray instru- 
ments on-board Chandra and XMM-Newton, while the X-ray non- 
detections are mostly due to ROSAT observations where the sources 
are located at large off-axis angles. For most of the sources with 



sensitive observations, X-ray information is available from pub- 
lished work (Gallagher et al. 2002, 2005; Xia et al. 2002; Maiolino 
et al. 2003; Ptak et al. 2003; Franceschini et al. 2003, hereafter F03; 
Satyapal, Sambruna & Dudik 2004; Ballo et al. 2004; Braito et al. 
2004; Imanishi & Terashima 2004; Iwasawa et al. 2005); for a cou- 
ple of sources, Chandra data have been retrieved from the archive 
and analysed using standard procedures. 

We find that AGN emission (with intrinsic X-ray luminosi- 
ties above fa 10 42 erg s" 1 ) is present in eight OHMs. Four of 
these OHMs are likely characterized by Compton-thick absorp- 
tion (UGC 05101, Maiolino et al. 2003; NGC 3690 in the galaxy 
pair Arp 299 system, Delia Ceca et al. 2002, Ballo et al. 2004; 
NGC 4418, Maiolino et al. 2003; Mrk 231, Gallagher et al. 2002, 
2005; Braito et al. 2004). However, Compton-thick absorption can- 
not be ruled out in some low signal-to-noise ratio sources and in 
Arp 220 (Iwasawa et al. 2001, 2005). 

To provide a more appropriate comparison with our Chandra 
targets, from the subsample of literature OHM with X-ray cover- 
age, we selected those within the same redshift interval, the same 
Loh range, or the same Lfir range as our seven OHMs observed 
by Chandra. The paucity of objects in the same redshift range 
(two, with only one X-ray detection) does not allow for any rea- 
sonable comparison with our sample; however, it highlights the 
relevance of the present work, where the redshift range of OHMs 
with X-ray constraints has been extended significantly. If we select 
the literature OHMs using the same interval of OH or FIR lumi- 
nosities as our targets, we find a larger number of objects (12 and 
nine, respectively). AGN emission is clearly present in three ob- 
jects (one Compton-thick, Mrk 231, and the other two Compton- 
thin, Mrk 273 and PKS B1345+145), corresponding to fa 25-33 
per cent of the chosen subsamples. LLAGN could still be present in 
some of the remaining sources, although the X-ray emission from 
these objects can be well explained by a thermal gaseous compo- 
nent from a starburst plus a hard, moderately absorbed power-law 
component due to the unresolved binary population (e.g., F03; Teng 
et al. 2005). 

As a further test, we compared the X-ray spectral results for 
our seven OHMs (see §4. 1) with those obtained by F03 with XMM- 
Newton for their five "bona-fide" lower redshift OHMs without a 
dominant AGN component (i.e., excluding Mrk 231). We found 
that our stacked spectrum is generally harder than the X-ray spectra 
of the F03 OHMs. This result can be ascribed to the different spatial 
resolutions of the Chandra ACIS and XMM-Newton EPIC cam- 
eras and different source-extraction regions. To test this possibility, 
we used a sample of three starburst galaxies observed by Chandra 
(Ptak et al. 2003) and XMM-Newton (F03). The soft thermal com- 
ponent appears spatially more extended than the hard component, 
i.e., it likely provides a stronger contribution to the broad-band 
X-ray spectrum if large source-extraction regions are chosen (as in 
the case of XMM-Newton spectra, extracted from circular regions 
of radius 20"). Motivated by these considerations, we rescaled the 
source-extraction regions of our OHMs to match those adopted by 
F03 (taking into account the different redshift range of the two sam- 
ples), obtaining a somewhat softer X-ray spectrum (r fa 1 .6 ± 0.3), 
consistent with F03 results. 

In the light of the results obtained from the study of the liter- 
ature OHMs and keeping in mind all the caveats described above, 
starburst emission with moderate absorption probably within the 
binary population is likely to explain the average X-ray properties 
of our seven OHMs (§4.1). 



5 SUMMARY 

We have analysed Chandra snapshot observations for a sample of 7 
galaxies selected from the Arecibo OHM survey at z ~ 0. 13 — 0.22 
and with FIR luminosities in excess of 10 11 L . The principal re- 
sults of this work are 

• All of the observed sources are X-ray weak; only one OHM, 
IRAS FSC 03521+0028 (z = 0.15) , was detected by Chandra, 
with 5 counts in the 0.5-8 keV band. 

• The X-ray emission is consistent with that from star formation 
(as also suggested by the comparison with most of the literature 
OHMs with X-ray constraints) and LLAGN; however, if an AGN is 
present, its contribution to the broad-band emission of OHM galax- 
ies is probably modest. 

• Under reasonable assumptions about the intrinsic X-ray spectral 
shape, the observed count distribution from stacking analysis sug- 
gests absorption of ~ 10 22 cm~ 2 . 
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